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Abstract 

We carry out a numerical and analytic analysis of the Yang-Lee zeros of 
the ID Blume-Capel model with periodic boundary conditions and its gener- 
alization on Feynman diagrams for which we include sums over all connected 
and non-connected rings for a given number of spins. In both cases, for a spe- 
cific range of the parameters, the zeros originally on the unit circle are shown 
to departure from it as we increase the temperature beyond some limit. The 
curve of zeros can bifurcate and become two disjoint arcs as in the 2D case. 
We also show that in the thermodynamic limit the zeros of both Blume-Capel 
models on the static (connected ring) and on the dynamical (Feynman di- 
agrams) lattice tend to overlap. In the special case of the ID Ising model 
on Feynman diagrams we can prove for arbitrary number of spins that the 
Yang-Lee zeros must be on the unit circle. The proof is based on a property 
of the zeros of Legendre Polynomials. 
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1 Introduction 



More than 50 years ago C.N. Yang and T.D. Lee started a theory |lj for studying 
phase transitions based on the partition function zeros . The same authors proved 
the celebrated circle theorem in |2j. It assures, in particular, that the zeros of 
the s = 1/2 Ising model partition function on a given lattice, occur only at pure 
imaginary values of the magnetic field. Defining u = e~^^'^ this implies that the 
zeros are on the unit circle in the complex u plane . The theorem is very robust 
and does not depend on details of the lattice like space dimensionality, number 
of nearest neighbors and topology. Later, the circle theorem was generalized to 
higher spin models without and with|3] extra symmetric potentials, as well as, to 
ferroeletric models and some Heisenberg models jHiinilZI- In all those models the 
spin variables are defined on a given static lattice and therefore they are the only 
relevant degrees of freedom. The proofs are based on the fact that the partition 
functions can be written as polynomials of a specific form 2J. Since addition of 
polynomials mix up their zeros in a rather non trivial way it is unexpected that the 
circle theorem will hold for spin models defined on dynamical random lattices which 
are themselves also degrees of freedom and must be summed over altogether with 
the spins. However, numerical results for the Ising model defined on 2D dynamical 
random lattices of planar [HI [HI HO] and torus topology [TT] indicate that the Yang- 
Lee zeros do lie on the unit circle although no circle theorem is known in the case 
of dynamical lattices. It is not even clear whether the new polynomials are of the 
form assumed in [2]. In particular, the form of the polynomials used in are not 
kept by arbitrary linear combinations. The dynamical lattices used in |H]-fI] were 
suggested in ^2] in order to mimic the effects of 2D fluctuating random surfaces in 
the continuum limit (2D gravity). The partition functions were generated by x 
Hermitean random matrices and the 2D gravity interpretation is achieved by fine 
tuning the coupling constants in the model and taking N ^ oo simultaneously. For 
the special case where the matrices become numbers (A^ = 1) although the 
2D gravity interpretation is lost, the models defined on the so called thin graphs, 
exhibit phase transitions with mean field exponents and strong similarities with 
models defined on the Bethe lattice ^Ij. In this case the partition function is 
obtained by summing over spin variables and connected and non-connected regular 
graphs (non-connected Feynman diagrams). Once again, numerical results valid for 
finite number of spins ^1] and analytic results valid in the thermodynamic limit 
HH] strongly support a circle theorem for this type of random lattices but no 
explicit proof for finite n has ever been given. It is the purpose of this work to 
produce such a proof at least in the simplified case of one dimensional models where 
the random graphs become connected and non-connected rings (see fig. 1). We 
have chosen the ferromagnetic Blume-Capel (BC) model ^Tj since it is the simplest 
model with "would be" Z2 symmetry such that on one hand their zeros ,in any 
dimension, are known to satisfy the circle theorem in a certain parameters range |3] 
and on the other hand they have been shown in 2D to exhibit non trivial topological 
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properties like bifurcation for other range of the parameters ^H|- Besides, the BC 
model includes the s = 1/2 and s = 1 Ising models. 



2 The ID Blume-Capel model on non-connected 
rings 

In the spin one Ising model there are three spin states on each site : Si = —1, 0, 1 . 
The Blume-Capel (BC) model is a generalization of the spin one Ising model where, 
besides the magnetic field H which works as a source for dipole interaction, we now 
introduce a coupling A functioning as a source for quadrupole interactions such that 
we still have an "would be" Z2 symmetry by exchanging H —H. The energy and 
the partition function of the Blume-Capel model are given by: 

n n 

E = -JY.s^s, - HJ2s^- xY.i^- sf) (1) 

(ij) 1=1 i=l 

= Y^e^ E(,,> s, + y:,^^ [hs, + A{i- sf) ] 
{s^} 

We only work with the ferromagnetic case J > 0. The sum is over nearest 

neighbors sites. For the one dimensional model we assume periodic boundary con- 
ditions Sj = Sn+j- In this case the lattice becomes a ring with n sites (vertices) and 
n bonds. Throughout this work we use the following definitions : 



K = J/kT 

h = H/kT 
A = X/kT 



c = e 



-K 



X = e 



u = (3) 

A 



Although originally suggested in a magnetic context ^7], the BC model corre- 
sponds to a special case of the BEG model [1^] and it can be used to describe phase 
separation driven by superfluidity in He^^-'-He^^-' mixtures. Its phase diagram in 2D 
is reach, including tricriticallity. At a; = (A ^ —00) it becomes the s = 1/2 Ising 
model while for a; = 1 (A = 0) it corresponds to the s = 1 Ising model. 

Each configuration with p spins down, Uq spins zero, and Uab bonds connecting 
spins Sa to Sb will contribute to the partition function a factor: 

g[(n-p-no)/i-p/i+noA] ^[K{n+++n — )-Kn+-\ j^^-j 

That factor will be needed in order to define the model on Feynman diagrams 
below. We start with the following generating function of Feynman diagrams in a 
zero dimensional field theory: 
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G{h,c,A,g) = „^ ^ , 

J d(t)+d(t).d(t)oe-2[^-^-*>M 



where the couphng constant g will play the role of a counting parameter while the 
3x3 symmetric matrix Mab-, a, 6 = +, — , is to be determined from the Boltzmann 
weights in (j2)) for /i = = D . The zero dimensional fields (f)± , 0o will represent the 
spin states Si = ±1,0 respectively. 
The expansion : 

oo 

Gih,c,g) = Y.9^2T , (6) 

n=0 

defines the coefficients 

^r = ^((^VUe-V- + e^0^)") (7) 

If we interpret e^cp"^ , e'^cj)"^ and e^0Q as interaction vertices representing sites 
with up, down and vanishing spins respectively, the Z^'^ have an interpretation as a 
diagrammatic expansion a la Feynman. In figure 1 we draw them for n = 1,2, 3. The 
Z^l'^ correspond to the sum of all diagrams with a total of n interaction vertices such 
that each vertex has two lines attached and all lines must be connected as in figure 1. 
Each vertex must be of course of type up, down or zero and will carry a factor e^, 
and respectively. Thus, we end up summing over all spin configurations altogether 
with all connected and non-connected rings with n vertices. The superscript "nc" 
stands for non-connected since the expansion includes non-connected diagrams like 
the second one in the second row and the second and third ones in the third row of 
fig.l. Each internal line of the diagram connecting a vertex of type a to a vertex of 
type h corresponds mathematically to a Gaussian integral < (pa't'b >= see ©■ 
Therefore, according to the so called Feynman rules^, each diagram with p minus 
and no vanishing spin contains a factor 

^[{n-p-no)h-ph+noD] (g) 
a<b 

In order to make more contact with the statistical model and match the above factor 
with (jU we require < (/>a(^fe > to be proportional to the Boltzmann weights in (j2)) 
calculated at /i = = A. That is. 



For a demonstration of the Feynman rules in a simple case see j2(J) 
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Figure 1: The Feynman diagrams corresponding to Z'^'^ ^ Z'^'^ and Z.^ 
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If we choose the overall constant k = cj [(1 — c^)(l — c)] we obtain a simple form 
for Mah such that the generating function ^ becomes 



G(c, /I, A, 



/ c/0+ c/0oe "^^ 



With 



1 



2 



1 + <P^ + f<Pl + 2c0+0_ - 2 (1 + c) (0+ + 0_) 00 



and 



/ 



(l + c)(l + c2) 



The Gaussian integrals in (jH)) can be easily calculated implying: 



G{c, h,A,g) 



1/2 



Where 



(9) 



(10) 



(11) 



(12) 



PM = {l-c'Y{l-c)-g{l-c'){5: + A)+g' (1 + c)(l + c^) + AS - g^x (13) 
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A = e^ + e ^ ; x = xc 
Therefore, from (jH)) we can write down 2^^^ as: 



(i-c^r(i-c) 



1/2- 



(14) 



(15) 



or equivalently from (|7j): 



5oe 



(>2_ + e^<Pl 



It is instructive to relate Z'^^ given in or (fTHIl with the usual partition functions 
of the one dimensional BC model on a ring Z^ given in (j2)). The Feynman rules 
and our definition for the overall constant k allow us to write the diagrammatic 
expansion: 



■-znc 



(l-c2)(l-c) 



^ Zn{diagram) 
(diagram) 

diagrams '^^ ' 



(17) 



Where the sum is over all diagrams with n vertices, connected and non-connected 
. The factor Sn is a positive real number called symmetry factor of the graph and 
it comes from the sum over all possible ways of connecting the vertices among each 
other to build the corresponding graph. The number of such possibilities we call 
An. According to our notation, see (|7j), and the Feynman rules Sn = (2"n!)/74„. In 
order to obtain An it is convenient to label the two lines of the i-th vertex as (oj, 6j). 
For example, for the case of connected rings, like the first diagrams on each row of 
fig. 1, there will be a total of 2n labels, if we start with the label ai it can be linked 
with any of the labels aj or bj such that j ^ l- Thus, we have 2n — 2 possibilities for 
the first connection. The yet non-connected label bi can be linked with the labels 
of any vertex but the first and the j*'^ ones. Thus, there are 2(n — 2) possibilities. If 
we keep going until no label is left unconnected we end up with An = (2n — 2)!! and 
S'n (connected ring) = (^2ri^2y.\ ~ "^^^ symmetry factor of the non-connected 



diagrams can be similarly obtained. As an example, for n = 3 (third row in fig. 1) 
the expansion (fTTj) becomes: 



^3 



6 



+ 



Z1Z2 



+ 



48 



;i-c2)(l-c). 

The general expansion for arbitrary n can be obtained as follows. First, it is well 
known in field theory that the connected diagrams can be singled out by taking the 
logarithm of the generating function. Using that Sn = 2n for connected diagrams 
we have: 
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\iaG{h,c,A,. 



OO -| 

y — 



(19) 



^_^2n [(l-c2)(l-c 

Where on the right handed side we have only partition functions of connected rings. 
Exponentiating (fTT?|l and comparing with the expansion (jHI) we have the generahza- 
tion of ()18|) for arbitrary number of vertices: 



(c, h) 




In particular, compare with p7|) . the symmetry factors of the non-connected dia- 
grams are produced automatically. The expression (PT|) shows that a generalization 
of the circle theorem valid for Zn to sums like Z'^'^ is non trivial. 



3 Yang-Lee Zeros 

3.1 ID BC model on a connected ring 

Before we start with the ID BC model defined on Feynman diagrams we first study 
the zeros of the static lattice (one ring) model. 

The partition function (j2]) is a polynomial on each of the variables c, x and u. 
In particular, Zn is a polynomial of degree In in the "fugacity" u. Concerning the 
Yang-Lee zeros on the complex -u-plane, it was proven in j3] that they must satisfy 
the circle theorem for < x < 2 independently of lattice details like topology, 
coordination number or space dimensionality. The proof, reproduced below, is based 
on the decomposition [2j of a spin one particle in a cluster of two spin one-half ones: 
Si = (cTji -|- crj2)/2, where (Jii;ai2 = ±1. After such decomposition we can rewrite, 
up to an overall constant, Zn as a partition function of a bigger system with 2n spin 
one-half particles: 



Zn OC Z; 



2n 



J2 

{a-ot} 

where a a include all an and with i 



h 2" 

<al3> a=l 

1, ■ ■ ■ , n. The sum J2 



<al3> 



(22) 

runs over all 



nearest neighbors (n.n.), including spins in the same cluster. The new coupling is: 



(23) 



(In 2 — A) within a cluster. 
K/A between n.n. of different clusters. 

The factor In 2/2 comes from the double degeneracy of the S* = state which can 
be represented inside a cluster by the weight Wi^aa.o-p) = exp |ln2 {aaO'13 — 1 
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Since the circle theorem holds for any ferromagnetic partition function of the form 
(|22|) we conclude that for A < ln2 (0 < x < 2), see j3], the Yang-Lee zeros of the 
BC model will lie on the unit circle for any space dimension. 

To the best we know, even though we will be dealing with the exactly solvable 
one dimensional model, there are no exact analytic results about the Yang-Lee zeros 
of the BC model for A > ln2. The best we can do is to present what we found 
numerically for finite number of spins and deduce approximate analytic results along 
the lines of the work of fIT] for the one dimensional Potts model. The starting point 
for the numerical calculations of the zeros is the exact expression for Zn which can 
be obtained , for instance, by inverting the expansion (fT^ using (O, 



Zn 



2n 



In 



1 -c^)^(l -c 



= (A'^ + A'! + A^) /c" 
where X± and Aq are solutions of the cubic equation: 

X^-X^A + x) + X \A{1 -c)x + 1- c^l - x{l - c^){l - c) 







(24) 
(25) 

(26) 



Incidentally, we notice that the inversion of formula is a diagrammatic 
method alternative to the transfer matrix approach. The quantities A-t/c; Aq/c 
are precisely the eigenvalues of the transfer matrix of the ID BC model, although 
we never really used it in our derivation of which came from factorizing Psig)- 
We remark that the way we found the exact solution (j^^ is a combinatorial one, 
where the FejTiman rules take care of the correct combinatorial factors. In this sense 
it is similar to the original solution of the one dimensional Ising model by E.Ising 



Proceeding further, for the purpose of generating the exact partition functions 
for large n it turned out to be faster to use the expression instead of plugging, 
for fixed c and x, the solution of in ()25|1 . For other higher spin models s > 1 
the computational advantage of against is even bigger. Given the exact 
expressions for the partition function, we have numerically obtained their zeros using 
the software Mathematica. Among other checks we have used the exact analytic 
results for the zeros of the ID Ising model on one ring to confirm the accuracy 
of the zeros. As a general feature we have found numerically that the complex 
u = zeros never converge for non-zero temperatures to the positive real axis in 
agreement with the fact that no phase transition for real magnetic fields appears in 
one dimensional spin models with short range interaction. Furthermore, the zeros 
lie on the unit circle for < x < 2 and arbitrary temperatures < c < 1, thus 
confirming the analytic result of For any given x > 2 we have always been able 
to find a temperature after which the zeros departure from the unit circle as shown 
in figures 2a-2d. 
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Fig. 2a 
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Fig. 2b 
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Figure 2: Yang-Lee zeros of the usual ID Blume-Capel model on the complex e 
plane for x = 3 and temperatures c = 0.435; 0.450; 0.500 and 0.520 respectively in 
figures 2a, 26, 2c and 2d. 



As X increases the minimal temperature for which that happens decreases. As a 
general rule, increasing x or c or even both of them the zeros move in a such way 
that after a given value for the product xc = e^~^ they start bifurcating into two 
branches: one inside and the other outside the unit circle. A bifurcation of the 
Yang-Lee zeros has been also observed in the two dimensional BC model but 
in that case the bifurcation starts at ^ = instead of ^ = vr. Our numerical results 
show that in the interval < x < 1 the zeros lie on the unit circle^ Although we 
have not been able to deduce a rigorous proof of this fact analytically, we have 
approximate analytic results based on the approach of [21]. The main idea is that 
the zeros of appear whenever at least two transfer matrix eigenvalues have the 
same absolute value which on its turn must be bigger than the absolute value of the 
remaining eigenvalue, e.g., |A+| = |A_| > |Ao|. Unfortunately, differently from the 
Q-state Potts model treated in |2I] which is similar to the Ising model {Q = 2), the 
secular equation (pUj) is not easily factorizable and their explicit solutions are not 
very useful to match absolute values. We have to stick to approximate solutions. In 
order to keep some similarity to the Ising model (x = 0), we solve (j^Bj) as a series 
expansion for low x. That gives: 

^The numerical results presented in |18| at a fixed temperature in the 2D case are in approximate 
agreement with the interval < i < 1, where of that reference corresponds to 
comment below figure 1 of that reference. 
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Ao 



1 - C2)2x 2Ac(l - Cfx^ 



1 + C2 



Vl-c^g± + 



l + c 



2\3 



C + 



2c^ 



where 



A 



9± 



± 



A 



X + 0{x 



1/2 



(27) 
(28) 



(29) 



The singularity at g'^ = 1 is the Yang-Lee edge singularity (A 



pe 

|2 



: A_) of the Ising 
*^ and (?-_ = e^*^/ p. Therefore, 
= r{l/p, -e,c,x)r{l/ p 9,c, x ) 
1 {A = cosh(/i) = VT^), 
which leads us to the unit circle \e'^\ = 1 since < cosh(/i) < 1, we have |A+| = 
|A_|. Therefore, our assumption will lead in fact to zeros on the unit circle if 
S± = |A±| — |Ao| > 0, or equivalently S± = 3fJe(/± — /o) > where fi = InAj. This 
can be checked from (jTfjl and (j2HI)- Using our assumption g± = e^*^ we get : 



model (x = 0). Since g+g- = 1 we can write g+ = 
we write |A+p = r{p,9,c,x)r{p,—9,c,x) and |A_ 
for some function r{a,b,c,d). Thus, assuming p 

■,h\ 



-In 



1 + cy{i + c) 

(1 - C)3 



In x + xc 



5 cos a 



l + c 



,2^2 



+ O X 



Z.2 



(30) 



Where h = ia, consequently cos 6* = cosa/vT^-c?. The first term on the right 
handed side of (fHUjl is always positive and the second one shows why the zeros go 
to the unit circle for < x < 1. The third term is also in agreement with our 
numerical calculations since it is more negative, thus pushing the zeros outside the 
unit circle, for angles close to a = tt. Adding up the three terms we have checked 
that 5-1- > whenever x < x* where the minimum value for = 0.780(2) occurs 
at c = 0.220(5) and a = n. The fact that x^i^ 7^ 1 is an artefact of the low x 
approximation. In practice our numerical calculation of the zeros show that they 
lie on the unit circle for < x < 1 with great accuracy. Although that seems to 
be a weak form (temperature dependent) of the circle theorem, one can slightly 
redefine the BC model, as in such that the three configurations where all spins 
are on the same state become the zero temperature degenerate ground states. This 
redefinition (x — > xc = x) makes x temperature independent, however the results 
derived so far for the zeros on the complex u = plane are still the same. 

After this partial understanding of the conditions under which the Yang-Lee zeros 
remain on the unit circle, one might ask what can we say about the parameters region 
for which the Lee- Yang circle theorem is not valid. We do not have a nonperturbative 
approach to this question but we have made low temperature expansions (LTE) 
analogous to (jHUjl . First of all, we found useful to implement the redefinition x — > 
xc = X before taking c ^ 0. We emphasize that such redefinition does not change 
the position of the zeros . Solving by a power series on c we have: 
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2(Ax — l)x / 9\ /X 

^ - - + (,-e-.)(J.->) ^ + ° ) (31) 

^± = + + (32) 



Defining = pe*", imposing that 3fJe(/+ — /o) = or 3fJe(/_ — /o) = we obtain 
respectively: 



X (x^ cos a — 1) / 2 

P = P+ = X + — — ^ — — ^ c +0(0 

1 + — 2x"^ cos a ^ 

1 (x^ cos a — 1) / 2\ 

X x(l + X* — 2x"' cosa) ^ ^ 

Notice that, due to the H —H symmetry which exchanges /+ and /_ we have 
p+ = Assuming |A+| = |Ao| we are led to p = p+, in this case we have 

calculated the difference: 



- /-) = 21n£ + 3(Y"s« 1) , + o (c') ^ F{x,c,a) . (33) 
X* + 1 — 2x^ cos a ^ ^ 

Whenever F > we have |A+| = |Ao| > |A_|. If we had assumed p = p_ in 
calculating 3fJe(/+ — /_) we would simply exchange the sign of the result, i.e., 
3?e(/+ — /_) = —F and consequently for F > we would have |A_| = |Ao| > |A+|. 
Concluding, for F > the zeros may come either from |A+| = |Ao| > |A+| with 
p = p+ or I A- 1 = |Ao| > |A+| with p = p_ while for F < they come from 
|A+| = I A- 1 > |Ao| with p = 1. Therefore, the form of the function F shows, at least 
for low temperatures, why the zeros departure from the unit circle as x increases 
beyond x = 1. In figure 3 we plot p+ and p_ for c = 0.1 and x = 1; 1.1; 1.2 altogether 
with the numerical results for the zeros for n = 100 spins. 

We see that the zeros inside (outside) the unit circle basically come form |A_| = 
|Ao| ( |A+| = |Ao| ) while the zeros on the unit circle originate from |A+| = |A_|. One 
can check that F > for the cases of figure 3b and 3c but this is not true for figure 
3a which explains the disagreement with the analytic curves. The large deviations 
which appear in the region close to the edges of the curves of zeros are consequences 
of the bad behavior of the LTE around the Yang-lee edge singularity. The same 
problem occurs in the LTE of the simpler case of the s = 1/2 ID Ising model. In 
passing, we notice that the Yang-Lee edge singularities can be obtained by imposing 
the degeneracy of the two largest eigenvalues of the transfer matrix as in [22] • We 
have used those singularities as a double check on the position of the zeros at the 
edges of the arcs. For the number of spins used in figures 3 the Yang-Lee singularities 
basically overlap with such zeros. In particular, in figure 3c there are four Yang-Lee 
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Figure 3: Overlap plot of p+ (dark solid) ,p_ (light solid) and Yang-Lee zeros (dots) 
of the usual ID Blume-Capel model at c = 0.1 and x = 1.0 (fig.3a) , 1.1 (fig.3b) and 
1.2 (fig. 3c). In fig. 3a the zeros lie on the unit circle. In all figures we used n = 100 
spins (200 zeros). 



edge singularities coinciding with the four edges of the two arcs. Finally, the reader 
may ask whether a triple degeneracy is possible, i.e., |A+| = |A_| = |Ao|. In practice 
we never found this type of zero in our numerical results. After an exact analysis 
of the secular equation pHjl we have found a necessary but not sufficient condition 
: X > ^/YT~c. 

3.2 ID BC model on Feynman diagrams 

We begin with the simplest case of the ID Ising model on Feynman diagrams (x = 0). 
Only in this special case we have Aq = and the polynomial p3( (7) becomes a second 
order polynomial such that the Z^'^, as given in p5|) . turn out to be proportional to 
Legendre polynomials (P„) as follows: 



1 



.n=0 



( cosh.{h) 



C2)(l 



n/2 



cosh(/i) 



(34) 



where ~g = g V^T^ / [(1 - - cff^. 

It is known that polynomials which are orthogonal in the interval [a, h] must 
have their zeros inside such interval excluding the borders. Therefore, the zeros of 
the Legendre polynomials VniUi) = ; z = 1, ■ ■ ■ , are such that — 1 < ?/j < 1. The 
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Yang-Lee zeros of the ID Ising model on Feynman diagrams can thus, be obtained 
from the zeros of the Legendre polynomials: cosh(/ij) = \J\ — yj which implies 
that the magnetic field is pure imaginary {hj = iaj). So the Yang-Lee zeros must 
be on the unit circle on the complex u = plane. This is so far the only exact 
proof of a circle theorem for statistical models defined on Feynman diagrams valid 
for finite number of spins. For two dimensions we found an analytic proof for the 
Ising model on Feynman diagrams with cubic vertices in ^T] but valid only in 
the thermodynamic limit and inside the convergence region of the low temperature 
expansion. 

In order to single out the effect of the Feynman diagrams on the zeros, it is 
interesting to compare the zeros that we have just found with the ones of the usual 
ID Ising model with periodic boundary conditions (connected ring). We first notice 
that each zero of the Legendre polynomial furnishes two Yang-Lee zeros ±a"'^ given 
by 

cosaf = Vl - Vj . (35) 

For the usual ID Ising model (x = 0) we have Aq = 0, see (jSHl), and Zn = (A" + 
X^)/c^. The Yang-Lee zeros come from X± = re^*^'-' where = {2k — l)7i/{2n) with 
k = 1,2, - ■ ■ ,n. From (1^ we have A+ + A_ = A = 2 cosh(/i) and A+A_ = (1 - c^). 
Therefore r = a/1 — and the zeros must be on the unit circle hj = iaj with ±aj 
given by 



cosaj = Vl — cos . (36) 

Comparing (j35|) with (j36p we see that, while the zeros of the model on the static 
lattice could be analytically determined, the sum over Feynman diagrams have led 
us to more complicated zeros, though still on the unit circle. Only numerical results 
can be obtained in general for zeros of Legendre polynomials. In the thermodynamic 
limit n ^ oo we have asymptotic formulae for the zeros of Legendre polynomials, 
see e.g. f2^, at leading order : = cos [{ik — l)n/ (4n + 2)] with k = 1, 2, ■ ■ ■ , n. 
Using this formula it is easy to check that aj — is positive (negative) for zeros 
close to the positive (negative) real axis and it approaches zero at a close to 7r/2 in 
agreement with figure 4. 

Besides, the difference Uj — a^ is of order 1/n^ in general which means that in the 
thermodynamic limit the zeros of the ID Ising model on Feynman diagrams approach 
those of the model on the static lattice (one ring). At n = 70 spins we found already 
no visible difference between the zeros. This seems to be a general feature of the ID 
spin models defined on Feynman diagrams as we will see later. It is remarkable that 
the zeros of the usual ID Ising model with periodic boundary conditions are directly 
related, in the thermodynamic limit, with the zeros of Legendre polynomials. 

Next we analyze the general case of the ID BC model on Feynman diagrams 
for arbitrary x. Unfortunately, we have not been able to rewrite Zl!f' in terms of 
orthogonal polynomials for a; 7^ 0. Thus, for finite number of spins we had to find the 
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Figure 4: Yang-Lee zeros for the Ising model on a connected ring (empty squares) 
and on Feynman diagrams (filled triangles) at c = 0.1 and n = 20 spins. 



Yang-Lee zeros numerically. It turned out faster to calculate the Gaussian integrals 
in the expression (fTI)|) than to use the formula (fTH|l . We found the zeros with great 
accuracy for diagrams with up to n = 120 spins, i. e., up to 240 zeros. We have 
checked for different temperatures, different values of x = and various number 
of spins that whenever the circle theorem holds for the usual ID BC model it will 
also hold for the corresponding model on Feynman diagrams. So, there must be a 
generalization of the circle theorem for Feynman diagrams which is valid for finite 
size diagrams. 

Furthermore, as we approach the thermodynamic limit we found in all cases 
that the Yang-Lee zeros on Feynman diagrams always tend to the zeros on the 
static lattice (one ring), see figure 5 where the zeros of the two models are already 
close for tt, = 40 spins. This has been confirmed by an argument based on the saddle 
point equations as follows. After a redefiniton (pa — > \/n(j)a in the numerator of ()lfj|l 
and calculating the denominator we obtain: 

n+3/2|'i _ /I _ I, 

= ) ,V.n r / d<P^d<p.dcPoe-' (37) 



2«n!(7r/2)3/2 7^ 

Where S = Sg=o — In ^ with V = e^cj)"^ + e~^(f)'^ + e^0o- The saddle point equations 
defined by daS = are given by: 



0+ - (1 + C)0o + C0_ = 
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Figure 5: Yang-Lee zeros for the Blume-Capel model on a connected ring (empty 
squares) and on Feynman diagrams (filled triangles) atc=0.1,x = l.l and n = 40 
spins. 



</)„ - (1 + c)(j)o + C(j)+ 
/0o-(l + c) (0+ + 0-) 



2e' 



(38) 



2e\ 



V 



It is easy to combine the above equations to show that at all saddle points (s.p.) solu- 
tions we have Sg=o = 1 and V satisfies a cubic equation such that (1 — c^)(l — c)V/2 = 
A satisfies precisely fEB|) . Thus, the six solutions of the saddle point equations are 
doubly degenerated since they produce only three different values for S at leading 
order. Expanding about the saddle points 4>a = (i>a + Ca/ \/n we have up to quadratic 
terms S = S + {CaCbdadbSls.p^/^n. Integrating over the fluctuations d(a and using 
Stirling's approximation we have at leading order: 



2{1 -c^)VT^ 



V2 



nnc 



c2)(l 



E 

s.p. 



A 



s.p. 



detdadbS\s.p. 



(39) 



where the factor two in the numerator comes from the double degeneracy of the 
saddle point solutions and A^.p, can only be A+, A_ or Aq. If we compare the symmetry 
factors in ()21|) we notice that the more non-connected diagrams like the first one 
between brackets drop out at n ^ oo. However, at this limit, see ()21|) . the connected 
diagram is not the only one to survive. If we remember the solution of the usual 
connected partition function (|25p and look at (j39|) we realize that the effect of the 
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other partition functions of the same order as Zn is to produce the denominators 
^detc^a^bS'ls.p.. 

The sum in ()39p extends over all relevant saddle points. In order to single out 
which solutions of (jHS|) are relevant, a careful analysis of the absolute values of 
would be necessary which depends in general on the parameters c, e^, and is rather 
complicated if we want an exact result. Instead, we discuss the case of the Ising 
model (x = 0) and assume that the general case will be similar as our numerical 
results for finite number of spins seem to indicate. For the Ising model we know 
that Z"^^ are proportional to Legendre polynomials for which asymptotic formulas 
are known, see e.g. [21] so we can check the result. If a; = we have Aq = 0. We 
are thus left with two doubly degenerated saddle point solutions. We assume that 
they all contribute whenever |A-|_| = |A_| which means A-t = \J\ — e^*^. In this 
case (jnH) becomes. 




The above expression reproduces exactly, including the overall constant asymptotic 
formulas for Legendre polynomials at leading order [24 , compare with (jS3)) which 
shows a posteriori that all four saddle point solutions ( doubly degenerated ) must 
contribute. The Yang-Lee zeros correspond to cos a = y/l — cos [{Ak — l)7r/ (4n — 2)] 
with /c = 0, ■ ■ ■ , n as we mentioned before. Since A-t are the solutions of the equation 



ipUl) with X = 0, i.e., A-t = cosh(/i) ± ycosh^ — 1. li h and c are real both A± 
are also real and A+ > A_ . In this case det dadbS < at the second solution, see 
second term in PU)). Thus, the steepest descent trajectory only passes by the first 
saddle point, the second point is not a true saddle point. Once again we reproduce 
asymptotic formulas for Vn{y) for ?/ > 1. In particular, the free energy of the model 
on Feynman diagrams is the same one of the traditional ID Ising model on a ring 
up to the addition of a constant. 

As a final comment concerning the ID model Ising on Feynman diagrams, we 
mention that the condition det dad^S = at the saddle points gives rise presumably 
|15j to continuous phase transitions. In our case this implies that the Yang-Lee edge 
singularity of the Ising model defined on Feynman diagrams coincides with the same 
singularity in the usual ID Ising model, i.e., A-t = A_. 

For the general case of the ID BC model on Feynman diagrams, once again, the 
zeros can only appear whenever the absolute value of, at least, two eigenvalues Aq 
degenerate, which is the same basic equation for the zeros that we have obtained 
when discussing the traditional ID BC model on a ring. This equation gives p = p{a) 
in the complex = pe*" plane. The precise values for a will be fixed by the 



relative phases of the saddle points which are present in the factors J det dadbS\ 



\s.p.- 
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According to our numerical results for finite number of spins we believe that the 
will be close to their corresponding counterparts in the connected (traditional) ID 
BC model. 

3.3 Conclusion 

Concerning the usual ID Blume-Capel model defined on a connected ring (periodic 
boundary conditions) our numerical results have shown that their Yang-Lee zeros 
remain on the unit circle as far as either A < In 2 or A < . The first condition has 
been known for many years j3] to be sufficient, in arbitrary dimensions, to rewrite 
the Blume-Capel model in a such way that the original Lee- Yang circle theorem 
applies. The second one is new and it is not clear whether it only works in the one 
dimensional case. We have presented a perturbative proof of this second condition 
valid for small e^~^ . Furthermore, we have observed that as a general feature as we 
increase A or the temperature or even both of them, the zeros which lie originally 
on the unit circle tend to bifurcate as in the two dimensional case ^H] , although no 
phase transition is present for real magnetic fields and non-vanishing temperature 
in one dimension. At some point (A*, K*) two disjoint arcs appear. We have shown 
by a low temperature expansion that the existence of two curves can be understood 
from the degeneracy of the absolute value of the transfer matrix eigenvalues. One 
of the curves comes from |A+| = |Ao| while the other one is due to |A_| = |Ao|. For 
the zeros which remain on the unit circle the degeneracy condition is: |A+| = |A_|. 

Regarding the ID Blume-Capel model defined on Feynman diagrams, which in- 
cludes a sum over connected and non-connected rings, although its partition function 
is rather complicated, see (PT|) . its zeros seem to follow, specially for large rings, the 
same pattern of the corresponding model defined on a connected ring. In particular, 
the two previously mentioned conditions are apparently satisfied. We have not been 
able to furnish a rigorous proof of this fact in general but an argument based on 
the saddle point solutions show that the conditions for Yang-Lee zeros are the same 
degeneracy conditions of the last paragraph. We believe that this is a general feature 
of models defined on the type of dynamical lattice used here. This indicates that a 
generalization of the circle theorem for dynamical lattices might exist. In fact, in 
the special case of the ID Ising model A — > —oo we have proved that the Yang-Lee 
zeros of the model defined on Feynman diagrams do lie on the unit circle. This is 
the only rigorous proof we know valid for dynamical lattices and finite number of 
degrees of freedom. It is not clear whether the polynomials can be written in the 
form used in the original proof of the circle theorem in or whether this proof is 
a independent one. Remarkably, it turns out that the zeros in the variable e'^ of 
the partition function of the usual ID Ising model with periodic boundary condi- 
tions are directly related in the thermodynamic limit with the zeros of the Legendre 
polynomials. 

Finally, we remark that throughout this paper we have obtained the partition 
function of usual ID models with periodic boundary conditions (one ring) by taking 
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the logarithm of the generating function of all diagrams including non-connected 
ones. This well known result in diagrammatic expansions furnishes an alternative 
to the usual abstract but simpler transfer matrix approach in solving exactly ID 
models. It has already proven to be useful from the computational point of view 
since it became faster than using solutions of the secular equation and plugging back 
in the partition function. 
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